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A numerical technique for calculating effective actions of electromagnetic backgrounds 
is proposed, which is based on the string-inspired worldline formalism. As examples, we 
consider scalar electrodynamics in three and four dimensions to one-loop order. Beyond 
the constant-magnetic-field case, we analyze a step-function-like magnetic field exhibit¬ 
ing a nonlocal and nonperturbative phenomenon: “magnetic-field diffusion”. Finally, 
generalizations to fermionic loops and systems at finite temperature are discussed. 


1. Introduction 

The computation of effective energies and effective actions of a quantum system 
is a general and frequently occurring problem in quantum field theory. In a given 
classical background, virtual quantum processes such as vacuum polarization modify 
the properties of the vacuum and are responsible for new nonlinear and nonlocal 
phenomena. These phenomena can often be described by an effective action that 
takes the virtual quantum pr^esses into account. A prominent example is given by 
the Heisenberg-Euler actiontla, describing the nonlinear corrections to the Maxwell 
action which are induced by virtual loops of electrons and positrons. 

The Heisenberg-Euler action is the one-loop QED effective action, 

r~ = lndet(-i^ + 4-im), = - lndet(-(a + iA)^ + m^), (1) 

(here in Euclidean formulation), evaluated for a constant electromagnetic back¬ 
ground field, and thereby represents a valid approximation of low-energy QED for 
all backgrounds that vary slowly with respect to the Compton wavelength 1/m. 

Technical difhculties increase markedly if one wantsTo go beyond this approxi¬ 
mation. The standard means is the derivative expansionS; its application is however 
limited to backgrounds whose spacetime variation does not exceed the scale of the 
electron mass or the scale of the field strength. Moreover, the enormous prolifera¬ 
tion of terms restricts the actual computations to low orders. Only for very specif 
(and very rare) field configurations can the derivative expansion be summed upu. 
Numerical methods developed so far aim at the integration of the corresponding dif¬ 
ferential equation of the operators in Eq. (^); therefore, they have been applied only 
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to highly symmetric background fields up to nowB. Of course, also the brute force 
method of spacetime discretization (lattice) can be employed, which has its own 
shortcomings (finite-size problems, fermion doubling, problems with many scales). 

In this note, we advocate a recently developed numerical techniqueQ which is 
baseion the string-inspired worldline formalismB applied to QED with background 
fieldsO. The idea of the approach consists of rewriting the functional determinants 
contained in Eq. (|^ in terms of 1-dim. path integrals, which finally can be evaluated 
with Monte-Carlo techniques. In contrast to the techniques mentioned above, our 
worldline numerics makes no reference to the particular properties of the background 
field and is therefore applicable to a wide class of configurations. In the following, 
we describe the approach in more detail and list a number of examples of a basic 
nature. 


2. Basics of worldline numerics 


In the following, we first concentrate on scalar QED, starting with the unrenor¬ 
malized, but regularized Euclidean one-loop effective action in D dimensions in 
worldline representation^, corresponding to the second equation of ([Ij), 


rr£r'"[A]= [ 

1/A2 


dT 

— I 

T 


Af J Vx{T)e 

x{T)—x{0) 


-f dr^^+ie x-A(x(T))'j 


( 2 ) 


Here we encounter a path integral over closed loops in spacetime. Note that there are 
no other constraints to the loops except differentiability and closure; in particular, 
they can be arbitrarily self-intersecting and knotty. Introducing the Wilson loop 


—ie r dr x(t)-A(x(t)) . r , 4/ x 

W[A{x)]=e 0 


Eq. (H) can be rewritten in a compact form: 


T-'Scalar 
^ eff 
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dT 


(47r)^/2 


Xo 


T{D/2) + 1 


\W[A]), 


1/A2 


( 3 ) 

( 4 ) 


Here we also have split off the integral over the loop centers of mass, Xq := 
(1/T) dTx^{T), and ((... ))a; denotes the expectation value of ( ... ) evaluated 
over an ensemble of x loops; the loops are centered upon a common average po¬ 
sition Xq (“center of mass”) and are distributed according to the Gaussian weight 

exp[-/o^rfr x]- 

Crucial for the numerical realization is the observation that substituting t =: Tt 
and introducing unit loops y, 

y{i) ■= tG[0,l], 


( 5 ) 
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renders the the Gaussian weight independent of the propertime: i^(T)/4 = 

fgdty'^(t)/4 (the dot “ ' ” denotes always a derivative with respect to the argu¬ 
ment). Now, the expectation value of W[A] can be evaluated over the unit-loop 
ensemble y, 


{W[Aix)]), = {W[VTAixo + VTy)])y, ( 6 ) 

where the exterior T-propertime dependence occurs only as a scaling factor of the 
gauge field and the unit loops in its argument. In other words, while approximating 
the loop path integral by a finite ensemble of loops, it suffices to have one single 
unit-loop ensemble at our disposal; we do not have to generate a new loop ensemble 
whenever we go over to a new value of T. 


2.1. Renormalization 


The effective action is renormalized by adding counterterms to Tea in order to absorb 
the strongly cutoff-dependent parts. The counterterms are determined by a set of 
physical constraints such as the vanishing of Tea for a vanishing background and the 
value of the gauge coupling in soft processes at a certain scale. To be specific, we 
use standard Coleman-Weinberg renormalization conditions for the complete action 
r, containing the bare Maxwell action and the one-loop contribution of Eq. (Q), 


r[F = 0] = 0, 


dr 

S(W) 


fv4=mV2 


1 

4^’ 


( 7 ) 


where we identify y with the scale of soft photons measuring the Thomson cross 
section, i.e., y/m 0. In the examples of Sect. 3, we will only encounter the field 
strength F which is renormalization group invariant, since it is scaled by the cou¬ 
pling. The substitution F e^iF^ reexpresses the later results in terms of physical 
coupling and field strength, where cr and Tr are the renormalized quantities at the 
scale y mentioned above. 

From a technical viewpoint, the strongly A-dependent terms can be isolated us¬ 
ing a heat-kernel expansion of the Wilson loop (which can be performed for arbitrary 
backgrounds), e.g.. 


(w[VtA{xo + VTy)]^ = 1 - Lt^ Fy4A]{xo) Fy4A]{xo) + 0{T^) . ( 8 ) 

The counterterms required are in one-to-one correspondence to the terms of the 
small-T expansion, encoding the short-distance physics; the number of necessary 
counterterms depends on the dimension of spacetime. While this renormalization 
program is well under control analytically, the numerical renormalization is compli¬ 
cated by a further problem: evaluating {W) with the aid of the loop ensemble does 
not produce the small-T behavior of Eq. (H) exactly, but, of course, only within the 
numerical accuracy. Unfortunately, even the smallest deviation from Eq. (|^) will 
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lead to huge errors of the effective action, because it induces an artificial singular 
behavior of the propertime integrand. 

Our solution to this problem is to fit the numerical result for (W) to a poly¬ 
nomial in T in the vicinity of T = 0, employing Eq. as a constraint for the 
first coefficients. In other words, we insert the analytical information about the 
short distance behavior of {W) for small T explicitly. This fit not only isolates 
the strongly cutoff-dependent parts which subsequently are subject to the standard 
renormalization procedure, but also facilitates a more precise estimate of the error 
bars. Finally, employing this fitting procedure only close to T = 0, the infrared 
behavior {T oo) oi the integrand remains untouched, and our approach is imme¬ 
diately applicable, also in the case m = 0. 

2.2. Numerical simulation 

The numerical computation of the effective action can be summarized by the fol¬ 
lowing recipe: 

(1) generate a unit loop ensemble distributed according to the weight exp[— dty^/4], 
e.g., employing the technique of normal (Gaussian) deviates; 

(2) compute the integrand for arbitrary values of T (and xq)\ this involves the 
evaluation of the Wilson loop expectation value for a given background; 

(3) perform the renormalization procedure; 

(4) integrate over the propertime T in order to obtain the Lagrangian, and also 
over Xq for the action. 

There are two sources of error which are introduced by reducing the degrees of 
freedom from an infinite to a finite amount: first, the loop path integral has to be 
approximated by a finite number of loops; second, the propertime t of each loop has 
to be discretized. Contrary to this, the spacetime does not require discretization, 
i.e., the loop ensemble is generated in the continuum. 

Approximating {W) of Eq. (^ by an average over a finite number of loops, 
the standard deviation provides an estimate of the statistical error. Approximating 
the loops by the finite number ni of space points results in a systematic error that 
can be estimated by repeating the calculation for several values n/. The number 
ni should be chosen large enough to reduce this systematic error to well below the 
statistical one. It will turn out that the choice = 1000 and ni = 100 for D = 3 
{ni = 200 for 0 = 4) yields results, for the applications below, which are accurate 
at the per cent level. 

2.3. Loop Clouds 

Beyond any numerical efficiency, worldline numerics also offers an intuitive approach 
to effective Lagrangians or functional determinants: Formula (^ together with (^ 
provides for a descriptive interpretation of the quantum processes. For a given 
background, the effective-action density receives contributions from all values of the 
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Figure 1: (a/b) small/large-propertime contribution to the effective-action density, 
picking up small/large-scale information of the background field (UV/IR). 


propertime. For instance, for small propertimes, the size of the ensemble of unit 
loops is also small by virtue of the propertime scaling in Eq. (D> VTy. Hence, this 
loop cloud picks up small-scale information about the background field. By contrast, 
for large propertimes, the loop cloud becomes bloated and receives information 
about the behavior of the background field over large scales. This zooming in or 
out of the quantum vacuum is illustrated in Fig. ^ 

Analyzing vacuum polarization in a certain background with this concept of 
loop clouds thereby visualizes the nonlinearities and nonlocalities of the quantum 
effective action in a vivid way. 

3. Examples 

3.1. Constant magnetic background field 

Let us first investigate the efficiency of our numerical loop approach to the scalar 
functional determinant for a constant magnetic background field B. For this case, 
the Wilson loop expectation value is exactly known and independent of D: 

{W[A]) = for H = const. (9) 

Integration over the propertime leads us to the effective-action density, which we 
plotted in Fig. (||) for I? = 3 and D = A. The agreement between the exact and 
the numerical results is satisfactory; the exact results lie well within the error bars, 
which contain the statistical errors of the Monte-Carlo calculation as well as those 
of the fitting procedure near T = 0 (the latter actually improve the error estimate, 
since our exact knowledge about the T = 0 behavior has been inserted). Our 
approach is able to cover a wide range of parameter values: in H = 3, even the 
massless limit can be taken without problems. In D = A, the typical logarithmic 
increase in the strong-field limit with a prefactor proportional to the P function is 
visible. 

3.2. Step-like magnetic field: quantum diffusion 

As a specific example for a field configuration that cannot be treated in a derivative 
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Figure 2: One-loop effective-action density 4ff in units of (-B/(47r))^/^ in Z? = 3 (a) 
and D = 4 (b) versus m?/B for the case of a constant magnetic background field. 
The analytically known exact results (solid lines) are compared with the numerical 
findings (circles with error bars). 


expansion but inherently involves nonlocal aspects, we consider a time-like constant 
background field B, resembling a step function in space, 

Bix,y) = -e{x)B, A{x,y) = e{x) ^{y,-x) B , (10) 

working for simplicity in H = 3. In our approach, discontinuities do not induce 
(artificial) singularities, but are smoothly controlled by the properties of the loop 
ensemble. This loop cloud has finite extension and slowly varying density; while 
running with its center of mass towards and across the step, that part of the volume 
of the loop cloud which “feels” the magnetic field increases smoothly. 

In Fig. H, we plot our numerical result for the effective-action density across 
the step (for details, seeO). As expected, the effective-action density is nonzero 
even in the region x < 0 where the background field B{x) vanishes. In order to 
characterize this quantum diffusion of the magnetic field further, we fit the effective- 
action density CeS in the zero-field region using the ansatz“ 

leS = ^ exp(-/3i m |x| - 132 Vb |x|). (11) 

In fact, the numbers ~ 3.255 and P 2 — 0.7627 fit the data with a high accuracy. 
The occurrence of two different diffusion lengths ~ 1/(/3 ito) and ~ 1/{(32'/B) may 
come as a surprise, since one might have expected that the only scale in the field-free 
region is given by the mass m. However, since the loop cloud is an extended object, 
it will “feel” the magnetic field even at a large distance; hence, the magnetic field 
is also a valid scale in the field-free region. This is precisely a consequence of the 
inherent nonlocality of the quantum processes. Moreover, the diffusion phenomenon 
appears to be also nonperturbative in the same sense as the Schwinger mechanism 


•^The spacetime coordinates occuring here are, of course, identical to the center-of-mass coordinates 
xq of the worldlines (cf. Sect. 2); we dropped the subscript for simplicity. 
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Figure 3: (a) Diffusion profile of the effective-action density in the vicinity of the 
magnetic step [x = 0) in units of {B /for various mass-to-field-strength 
ratios, (b) A logarithmic plot reveals the exponential nature of the diffusion depth. 
For a check of the propertime-continuum limit, three different loop ensembles have 
been used: ni = 100 for set A, m = 75 for set B, m = 50 for set C. 


of pair productiorc. This is suggested by the functional form Eq. (|Tl|), which cannot 
be expanded in terms of the coupling constant, being rescaled in the field (only an 
expansion in terms of the square root of the coupling constant is possible). 


3.3. Constant magnetic field at finite temperature 

The presence of a heat bath can be taken into account within the worldline frame¬ 
work by using the imaginary-time formalism: the Euclide an tim e direction is com- 
pactified to a circle with circumference /3 = l/temperaturell3tdl.^ This changes the 
boundary conditions of the worldline; the loops do not have to close trivially in 
spacetime, but can now wind around the spacetime cylinder n times: 



x{T)—x{Q) 


00 

E 

n= —00 


Vx. 


x{T)=x(0) ,xo{T)=xo{0)+nl} 


( 12 ) 


Choosing a gauge where the time-like component of the gauge field is time inde¬ 
pendent, a Poisson resummation of the n-sum disentangles the winding around the 
cylinder from the trivially closed loops for any background, and the remaining path 
integral runs over nonwinding loops only. The numerical problem thereby reduces 
to the zero-temperature case with an additional summation over Matsubara modes, 
which can be performed to a high numerical accuracy. 

Confining ourselves to the constant-magnetic field case, tl^^act result for the 
finite-temperature contribution to the Wilson loop EV reads Ej’Ed 


(W[A])f = 2 


TB 


^ sinh Ti? 

n—1 


(13) 


^The propertime T should not be confused with the temperature, being characterized by its inverse 
0 in the following. 
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Concentrating on D = 3 scalar QED, we compare the exact results with the numer- 
ical estimates in Fig. The propertime integrand (plotted for m = 1//3, = B) 

as well as the effective-action density (plotted for m = 0) can satisfactorily be 
reproduced by the numerical approach. The drop-off of the integrand for small 
propertimes in Fig. ^(a) can be understood in terms of the loop clouds: for small 
T the loop cloud is simply too small to “see” the compactness of the Euclidean 
time direction. Furthermore, the Stefan-Boltzmann law for the free energy (= — 
effective-action density) ^ C(3)//3^ is rediscovered in the B = 0 limit. 

(a) propertime integrand (b) 




Figure 4: (a) Propertime integrand of the finite-temperature contribution to the 
effective-action density in D = 3 scalar QED (m = 1//3, = B). (b) Effective- 

action density contribution at finite temperature in units of 2/(47r/3^)^/^ versus 5/3^ 
in the massless limit. 


3.4. Fermions in the loop 


The worldline formalism allows for an elegant inclusion of fermions by represent¬ 
ing the Dirac algebraic elements by path integrals over closed loops in Grassmann 
space. Unfortunately, this cannot be translated into a numerical algorithm, since 
the kinetic action for the Grassmann loops is not positive definite; this prohibits a 
Monte-Carlo algorithm based on importance sampling. Instead, we choose to main¬ 
tain the Dirac algebraic formulation, which leads to the following representation of 
the spinor part of Eq. (H) : 


^spinor 1 1 [ ^ 

2 (47r)^/2 J T ' 

1/A2 

where the “spinorial” Wilson loop is given by 


-m^T 




IU,pi„[A] = IU[A] X TrPr exp 




(14) 


(15) 


Obviously, the spinorial Wilson loop factorizes into the usual one and a trace of a 
path-ordered exponential of the spin-field coupling. Although this path ordering is 
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difficult to handle in analytical computations, our numerical algorithm can easily 
deal with; this is because the path(=loop) is discretized and each point is visited in 
a path-ordered manner anyway. As a further simplification, the exponential of the 
spin-field GQupling can be projected onto an orthogonal basis of the Dirac algebra 
(see, e.g., H). 

Concentrating again on the constant-magnetic-field case, the analytically known 
result for the spinorial Wilson loop is 

(f^spi„[A]) = BTcothBT = coshST, (16) 

where the cosh term arises from the trace of the path-ordered exponential in Eq. ([l^ . 
Comparing the exact result with the numerical one in Z? = 3 as depicted in Fig.|^, 
the agreement in the parameter range B < is again satisfactory. Only in the 
strong-held limit does the algorithm become unstable owing to technical reasons: 
here the exponential increases of the cosh and the sinh function do not cancel nu¬ 
merically as precisely as they should. This technical problem can be solved by Hrst 
matching the exponential increases of the denominator and the numerator to each 
other and then performing the division. Whether this matching can be numerically 
implemented in a background-held independent way has still to be investigated. 
Beyond this, there are no problems of any fundamental kind limiting the applica¬ 
tion of the present approach to fermion determinants, in contrast to no-go theorems 
faced by lattice formulations. 


^spinor 



Figure 5: Effective-action density of spinor QED in D = 3 for the constant- 
magnetic-held case (units as in Fig. ||). 
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4. Conclusion 

Our results of the numerical worldline approach to QED effective actions are en¬ 
couraging: the successful description of the background conditions selected here 
indicates that the method can be applied to a very general class of backgrounds 
without the necessity of assuming slow variation, smoothness, or symmetries. In 
particular, it should be stressed that the algorithm makes no reference to the prop¬ 
erties of the background field. 

Beyond being an efficient method of computation, our approach also offers a 
vivid picture of the quantum world: consider a spacetime point a: at a propertime 
T; then, the loop cloud is centered upon this point x with Gaussian “density” 
and “spread”. Increasing or lowering the propertime T corresponds to bloating or 
scaling down the loop cloud or, alternatively, zooming out of or into the microscopic 
world. The effective-action density at each point x finally receives contributions 
from every point of the loop cloud according to its Gaussian weight and averaged 
over the propertime. This gives rise to the inherent nonlocality and nonlinearity 
of the effective action, because every point x is influenced by the field of any other 
point in spacetime experienced by the loop cloud. 

We would like to conclude with some further remarks and caveats. Of course, 
the Monte-Carlo procedure requires that the loops are distributed in Euclidean 
spacetime in order to guarantee the positivity of the action. This can lead to 
difficulties (also shared by many other approaches) when the method is applied to 
the computation of Greens functions with Minkowskian momenta. However, this 
does not pose a problem for the inclusion of Minkowskian electric fields which can 
simply be taken into account by an imaginary Euclidean field strength. For instance, 
pair production in inhomogeneous electric fields can be studied in this way. Finally, 
the analogy between worldline path integrals and Brownian motion implies that 
the loops in our loop cloud should have Hausdorff dimension 2; on a computer, 
this can, of course, never be realized. However, we expect that this may lead to 
sizeable errors only in the extreme case of a massless or very light particle in the 
loop and a background field with fluctuations on all scales; only then are the large¬ 
propertime contributions not suppressed, and the coarse structure of the large loop 
clouds are insensitive to small-scale fluctuations. This failure will manifest itself in 
the practical impossibility of performing the propertime-continuum limit. On the 
other hand, if the continuum limit can be taken, the loops will be sufficiently close 
to Hausdorff dimension 2. 
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